plot (1 - (R^2/2 + R + 1)*e^(-R)) R from 0 to 10